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The t-J Hamiltonian is one of the cornerstones in the theoretical study of strongly correlated 
copper-oxide based materials. Using the density matrix renormalization group method we calculate 
the phase diagram of the one-dimensional t-J chain in the presence of a periodic hopping modulation, 
as a prototype of coupled-segment models. While in the uniform ID t-J model near half-filling 
superconducting (SC) state dominates only at unphysically large values of the exchange coupling 
constant J/t > 3, we show that a small hopping and exchange modulation very strongly reduces the 
critical coupling to be as low as J/t 1/3 - well within the physical regime. The phase diagram as 
a function of the electron filling also exhibits metallic, insulating line phases and regions of phase 
separation. We suggest that a SC state is easily stabilized if t-J segments creating local spin-singlet 
pairing are coupled to each other -- another example is ladder system. 

PACS numbers: 74.20.-z,71.lO.Fd, 71.45.Lr,74.25.Dw 


Introduction — Since its introduction at the end of the 
1980’s, the t-J model Hamiltonian [I| has formed one of 
the cornerstones in the theoretical study of high tempera¬ 
ture superconductors (HTSs). It is a minimal low-energy 
model for the electronic and magnetic structure of the 
copper-oxide planes in HTSs and can be derived from the 
Hubbard model in the stroM coupling limit using second 
order perturbation theory [ 3 . Although the initial inter¬ 
est in the t-J model focussed on its two-dimensional (2D) 
realization, its one-dimensional (ID) version, which is of 
direct relevance to for instance doped spin-chain materi¬ 
als, provides a number of interesting phases that are also 
observed in the 2D context, such as a ™in gap phase 
and the occurrence of phase separationjljQ. The ID t-J 
model actually also displays a superconducting instabil¬ 
ity, even if all this is in a parameter range where the ratio 
of the exchange constant J and hopping t, near half-hlling 
is J/t ^ 3, which is not of relevance to real materials - 
the HTSs are rather in the regime where Jjt ^ 1/3. 
In the latter regime quasi-lD swtems e.g., t-J ladders 
can support superconductivity ptSI which is related to 
the substantial binding energy for two holes in even-leg 
ladders, giving rise to the presence of preformed Cooper 
pairs on the rungs. However, in a ID system this binding 
energy vanishes for physically relevant values of J/t ra¬ 
in this Letter we consider the ID t-J Hamiltonian in 
the presence of a periodic local modulation of t and J: 
within segments of length Si the hopping and exchange 
are constant, but the coupling between these segments is 
somewhat weaker, see Fig.[TJ In a solid such a local mod¬ 
ulation might for instance result from a periodic struc¬ 
tural modification or from electronic self-organization. 
Using the density matrix renormalization group (DMRG) 
method we establish that for hole-doped systems close 
to half hlling in the presence of a weak modulation not 
only a spin-gap forms but also that the holes pair-up al¬ 
ready for moderate values of J/t. A calculation of the 
Luttinger parameter shows the formation of a supercon¬ 
ducting (SC) state in a large region of the phase dia- 
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FIG. 1: (Color online) Schematic of the model for coupled t—J 
segments. Here t and J are the hopping integral and exchange 
interaction respectively with corresponding reduction St and 
SJ only for the bonds connecting end sites of the segments. 


gram, also in the physically relevant low-doping regime 
with J/t ~ 1/3. The transition point weakly depends 
on segment length and apart from superconductivity 
also metallic, insulating and phase separated regions are 
present in the calculated phase diagram. 

Model and Method — To investigate the effect of cou¬ 
pling of segments forming a ID chain, we consider the 
model Hamiltonian H = Hq-\-Hs where Hq is the uniform 
ID t-J Hamiltonian Hq = -tj^ia + H.c.^ -b 

^E,(s i ■ Si_|_i — where the sum is over all in¬ 

tegers i labeling the sites that make up the chain and Hs 
describes the local modulation and 

Hs = StY.{^ls,.aCj.S, + l,a + H.C.'^ 

3<r 

- dJ ■ Sj,S, + i- - -), (1) 

j 

where the sum is over all integers j. Hs represents the 
reduction in hopping 6t and exchange interaction dJ on 
the bonds connecting end sites of segments of length Si. 
The modulated bond consist of the rightmost site of 
segment and next site to it, which are indexed by j.Si 
and {j.Si + 1) respectively - a schematic of the model is 
shown in Fig. [1] In the Hamiltonian cj^ is the electron 
creation operator at site i with spin ct, Si is the spin-i 
operator, and ni the electronic number operator. In order 
to stay within the perturbative framework in which the 
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t-J model is derived from the Hubbard Hamiltonian, we 
retain the direct relation between St and 6J from second 
order perturbation theory: as J = /U, where U is the 

onsite Coulomb repulsion, we have that 


J — SJ f ^ ~ ^ 


( 2 ) 


In units of t, we are thus left with only two energy scales: 
St/t and J/t. The parameter SJ/J is fixed by the rela¬ 
tion in Eq. [2]. Clearly for St = 0 the model reduces to 
the regular ID t-J Hamiltonian. The electron density is 
denoted by n. 

The quantities of interest e.g., Luttinger parame¬ 
ter, spin gap, binding energy are calculated by DMRG 
method prl. [l^ on a lattice with upto 288 sites and fi¬ 
nite size extrapolations are performed to obtain them 
in thermodynamic limit (see supplemental material). 
As we have studied the system with segment length 
Si = 2,4,6,8,12 we choose the system size L = 
48,96,144,192, 240, 288. This allows us to have the num¬ 
ber of electrons N = nL and the number of segments in 
the system to be even so that the ground state corre¬ 
sponds to total spin S'J = 0. All the results are obtained 
with upto 2000 basis states, 15 sweeps and typical dis¬ 
carded weight ~ 10“® leading to error in energy of the 
same order. 

Spin gap and pair binding — We start our discus¬ 
sion with the calculations for the spin gap Ag in the 
system with coupled segments. The usual t-J model 
(i.e., 5t = 0) in ID has a spin gap phase, as is 
found in the studies by exact diagonalization in small 
systems Q, variational methods [2^ RG analysisand 
more recently by DMRG method [J. But in all the stud¬ 
ies spin gap phase appears only at low electron density 
and at large exchange J/t > 2, a parameter regime that 
is barely relevant to real materials. We show below that 
the spin gap phase can appear at small J/t and for weak 
hole doping (i.e., n —^ 1) when a non-zero St is intro¬ 
duced. 

The singlet-triplet excitation energy is given by the 
energy difference Ag = E{N,Sf = 1) — E{N,Sf = 0). 
Here E{N, S'J) is the ground state energy with quantum 
numbers N and total ^-component of spins S'J. At half¬ 
filling, n = 1, the uniform t-J model corresponds to a 
Mott insulator with zero spin gap. But if St is switched 
on, a spin gap Ag forms and increases monotonically, see 
Fig. [2](a). The gap reaches its maximum in the limit 
St/t = 1, where it corresponds to the finite size gap of 
decoupled Heisenberg segments. However, the situation 
changes drastically when holes are doped into the sys¬ 
tem, see Fig. [2jb) and (c). The overall spin-gap is much 
reduced but kept to be finite, and reaches its maximum 
for small St/t and gradually vanishes for larger values of 
the hopping modulation. Fig. [2Kb) shows the spin gap 
As for J/t = 0.8 at different St in the limiting density 
n —>■ 1 which corresponds to putting only two holes in 
different system sizes and extrapolating Ag to thermo¬ 
dynamic limit (see supplemental material). The calcu- 




FIG. 2: (Color online) The spin gap Ag and binding energy 
Ab at J/t = 0.8 as a function of St for (a) half hlling n = 1, (b) 
in the limit of half filling n —>■ 1 i.e., 2 holes in the system and 
(c) n = 23/24. The empty circles (red), filled circles (green) 
and stared points (blue) connected with solid lines represent 
the spin gap for segment size Si = 2, A and 6 respectively. 
The corresponding circles and points connected with dotted 
lines represent binding energy. 


lations were done for segment sizes Si =2, 4 and 6. 
Surprisingly, a very small St = 0.01 is already enough to 
produce a sizable spin gap in the system. The spin gap 
increases first as in half-filled case, but then decreases to 
zero for larger St. It is interesting that the maximum of 
the spin gap is not really sensitive to the segment size. 
We notice that the occurrence of maximum of Ag shifts 
to higher St if we increase segment size Si. It is related 
to slower reduction of (electron or hole) bandwidth by 
St for larger Si. Also shown in the hgures is the pair 
binding energy, which measures the stability of pairing 
and is defined as Af, = E{N ± 2, = 0) -I- E{N, S'J = 

0) — 2E{N ± 1, 5'J = ±1/2). The plot for At is shown in 
Fig. |2Kb) for Si = 2, A and 6 with the same points and 
colors as Ag but connected by dotted lines. The nega¬ 
tive values of A^ indicates the stability of pairing and a 
relation Ag = |Ab| in a metallic (M) regime suggests an 
occurrence of spin-singlet SG state. The spin gap and 
pair binding energy at a finite hole density of n = 23/24 
are shown in Fig. me). Gompared to the limit of half fill¬ 
ing n —t 1, their magnitudes are smaller but the over-all 
behaviors are almost unchanged. 

Superconductivity and Phase diagram — The DMRG cal¬ 
culations allow us to map out the phase diagram for the 
modulated ID t-J Hamiltonian including also potential 
states. It should be noted that the uniform t-J model 
shows a transition from M to SG state at J/t > 2 in low 
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FIG. 3: (Color online) (a)-(c) Phase diagrams on n—Jjt plane at 5t/t = 0.2 for S; = 2, 4 and 6 systems respectively containing 
different phases: Superconducting (SC), Metalic (M) and Phase separated (PS) as indicated. The points (empty circles) 
represent M to SC transition and the point (star) indicates density n —>■ 1 (system with 2 holes). The vertical lines at different 
commensurate fillings represent insulating line-phases. 


electron density regime. An even higher J/t is necessary 
to get SC phase with increasing electron density 0. In 
order to examine the possibility of superconductivity in 
system with coupled segments at different electronic den¬ 
sity n, we determine the metal-SC transition by calculat¬ 
ing the Luttinger parameter: Kp < 1 for a paramagnetic 
metal and Kp > 1 for a superconductor. In the ther¬ 
modynamic limit the Luttinger parameter is determined 
from the slope of the structure factor for the density- 
density correlation at wave vector fc —?> 0 limit [l3, H] • To 
avoid difficulties in the calculation of real-space density- 
density correlations at large distances and the subsequent 
Fourier transform at small fc, we calculate structure fac¬ 
tor directly in momentum space [Tsj and extract Kp in 
thermodynamic limit by finite size scaling (see supple¬ 
mental material). This requires calculating Kp for differ¬ 
ent system sizes (upto 288 sites). 

We determined the crossing point Kp = 1 for fixed 
different densities n, segment lengths Si, and modula¬ 
tion strengths 6t, that separates the M and SC phases 
as a function of J/t. The critical value is denoted by 
J^/t hereafter. At very large J, the system separates 
into electron- and hole-rich region and the onset of phase 
separation (PS) has been determined by inverse com¬ 
pressibility (see supplemental material). The resulting 
n—J/t phase diagrams containing PS, M and SC phases 
are shown in Fig. [3Ka)-(c) for 5t/t = 0.2 and S'; = 2, 4 
and 6, respectively. The (blue) star symbol represents 
the critical point jt in n —> 1 limit. It is remarkable 
that close to half filling (even for a system with only two 
holes i.e., n —>■ 1) J^/t is significantly reduced for all 
the cases of three SiS as shown in Fig. [3l Apparently, 
the critical value can be as small as J^^/t ^ 0.35 for a 
relatively weak modulation 5t/t = 0.2. 

This observation motivates to study how the critical 
point K/t changes as a function of St/t. Since a large 
reduction of J‘^/t can be acquired near half filling, we fix 


the filling at n = 23/24 and estimate the critical point. 
The critical value jt at n = 23/24 is shown in Fig.|4]for 
Si = 2,4,6,8,12. For each segment length even a weak 
modulation bt causes a large drop in jt from ^ 3.2 
at the uniform case St = 0. Depending on the segment 
size the lowest J^/t is reached for 6t/t between ~ 0.1 and 
~ 0.3. The calculations imply that an increase of St helps 
electrons in each segment to form pairs which can move 
easily unless St/t becomes too large, which then tends to 
localize electrons on individual segments for large St/t. 
This explains the steep minimum and then slow rise in 
K/t. At small St/t, the transition point K/t increases 
with Si because the pair formation gets weakened as the 
number of sites in a segment becomes larger. Ultimately, 
at very large Si, this situation will correspond to the 
uniform t-J model. On the other hand at large St/t, 
the pair formation is most effective but the movement of 
pairs is restricted due to the narrowing of bandwidth and 
a high J°/t is given accordingly. An intriguing thing is 
that for Si = 12 the strong reduction of K/t is obtained 
despite K/t ~ 3.2 for any St/t in the A; —>■ oo limit. 
It means that a system with Si = 12 is still far from 
the uniform t-J model. We also notice a sharp decrease 
in J'^/t as we increase 5”; = 2 to 12 initially at large 
St/t and then almost saturates. To explain this behavior, 
we calculate the spin gap of decoupled segments (i.e., 
St/t = 1) of different sizes starting from Si = 2. As 
expected, we see a sharp decrease of spin gap with Si as 
shown in the inset of Fig. |4l Plots in Fig. [4] also explains 
the appearance of finite spin gap discussed in Fig. [2Dc) for 
density n = 23/24 at J/t = 0.8. The critical J‘^/t where 
the M to SC transition occurs is lower than J/t = 0.8 for 
some range of St/t and for different Si. This means the 
system develops a finite spin gap entering into SC phase. 
Insulating line phases — Unlike the uniform t-J model, 
we notice some insulating regions (vertical lines in Fig. [31) 
appearing at various commensurate Filings. This type of 
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insulating region at particular commensurate fillings has 
been discussed in ladder systems [l^. Due to the restric¬ 
tion in hopping as we introduce 5t, the electrons may 
form super structure in each segment which leads to in¬ 
sulating behavior of the system. The insulating to M 
transition at n = 1/2 which is shown in Fig. [21Ia) for 
Si = 2 system has been discussed in terms of melting of 
Mott insulating states[2l|. We briefly recap its intuitive 
understanding which will guide us to understand the in¬ 
sulating regions at other fillings and Si = 4, 6 systems. 

Let us first consider Si = 2 system with no exchange 
interaction i.e., J —>■ 0 in t/ —>■ oo limit. Then the tight 
binding model with 5t > 0 at n = 1/2 forms bonding 
and antibonding states with spin a i.e., t—dimer states: 
\t — dim)a = {\cr 0) ± |0 (j))/\/2 separated by a gap 
Ad = 2dt at each segment. Due to this t—dimerization 
one electron per segment is in the bonding state at this 
commensurate density of n = 1/2 and the system be¬ 
comes a Mott insulator. But if we increase J, a ten¬ 
dency to make singlet pairs is introduced at each segment 
due to J-dimerization which leads to a transition to M 
states. If we increase J further where J-dimerization 
dominates, the system with these singlet pairs becomes 
SC and finally goes to PS at large J. Such a melting 
of insulating states happen also in S'; =4 system at 
n = 1/4,3/4 and Si = 6 system at n = 1/6,3/6,5/6 
as shown in Fig. [Hljb) and (c), respectively. No effec¬ 
tive pairing with odd number of electrons in a segment 
at these electronic fillings causes the melting of insulat¬ 
ing states with increasing J. However, even number of 
electrons in a segment as in S; = 4 system at n = 1/2 
and S; = 6 system at n = 2/6,4/6 prevent such melting 
of insulating states by making effective pairing between 
themselves even with increasing J. So the insulating 
states continue to exist until the PS boundary at these 
hllings with increasing J as shown in Fig. [^b) and (c). 
These insulating line phases and metal-insulator transi¬ 
tions are captured by calculating charge gap defined as 
Ac = {E{N -h 2) -H E{N - 2) - 2E{N))/2 for different 
values of J. 

Conclusions'. We have shown using density matrix renor¬ 
malization group method that a metal to superconduct¬ 
ing transition can occur in a physically relevant regime of 
exchange interaction and hopping (i.e.,J/t ^ 1/3) in ID 
system with coupled t-J segments. A moderate modula¬ 
tion 5t/t in the hopping causes a superconducting state 
for electron densities close to half filling. As compared 
to uniform ID t-J model, this implies an order of mag¬ 
nitude reduction in critical exchange interaction to hop¬ 
ping ratio. We have presented full n^J phase diagrams at 
5t/t = 0.2 with different segment sizes containing besides 


phase separation, metallic and superconducting phases 
also insulating line phases at commensurate doping con¬ 
centrations. These results may be of relevance to certain 
quasi ID materials, including cuprates in which typical 
exchange coupling is J ^ t /^ Ig. An example is the 
spin-Peierls system CuGe 03 [l 7 | which consists of linear 
spin-1 Cu 02 chains with alternating exchange coupling 



FIG. 4: (Color online) This shows the transition from metallic 
to superconducting phase at density n = 23/24 as a function 
t 2 and for systems with segment sizes 2, 4, 6, 8,12 as indicated. 
In the inset we plot As for systems with decoupled segments 
of different sizes at J/t = 0.8. 


and hopping strength, corresponding to a S'; = 2 system 
in our language. The results above suggest that dop¬ 
ing with few holes may turn the system to be supercon¬ 
ducting. Such values of exchange coupling and hopping 
strengths can in principle also be achieved in ID ultracold 
fermionic quantum gases [l^ and with polar molecules [l^ 
on optical lattices. It will be most interesting to estab¬ 
lish whether in such quantum simulator experiments a 
superconducting state can be stabilized by a relatively 
benign modulation of the ID t-J Hamiltonian as we pro¬ 
pose here. In this Letter, the segments are open chain 
coupled linearly as a simplest ID case. However, any 
shape of segment can be allowed if it creates locally hnite 
spin excitation gap. The network between the segments 
would be also flexible. We think that this simple con¬ 
dition possibly makes a great help for superconducting 
material design. For example, the ladder t-J system is a 
special case of the coupled-segment systems. 
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Here we present a detailed plots of finite size scaling analysis to get spin gap, charge gap and 
Luttinger parameter in thermodynamic limit. Luttinger parameter has been used to determine the 
transition between metallic to superconducting phases whereas spin gap and binding energy indicate 
the stability of pairing of electrons. Insulating line phases are captured by finite charge gap and 
onset of phase separation is obtained from inverse compressibility. 


We use DMRG method[l|, 0 to study the system consisting of coupled t—J segments in presence of periodic 
hopping modulations on a lattice with total upto 288 sites and different segment size S';. We calculate each of the 
quantities of interest e.g., Luttinger parameter, spin gap, binding energy, charge gap in thermodynamic limit by finite 
size extrapolation. As we have studied the system with segment length S; = 2,4,6,8,12 we choose the system size 
L = 48, 96,144,192, 240, 288 and more in some calculations. This allows us to have the number of electrons N = nL 
and the number of segments in the system to be even so that the ground state corresponds to total spin S'J = 0. All 
the results are obtained with at least 500 basis states, 15 sweeps and typical discarded weight ~ 10“® leading to error 
in energy of the same order. 

We make use of Luttinger parameter Kp to find the metallic {Kp < 1) to superconducting {Kp > 1) transition as 
a function of exchange coupling J and electronic density n. This enables us to map out the n — J phase diagrams 
presented for Si = 2,4,6 systems. The spin gap and binding energy indicate the stability of electron pairing in 
superconducting state. Insulating line phases at different commensurate hllings are captured by charge gap calculation. 
Finally inverse compressibility has been used to estimate the phase separation boundary. 


SPIN GAP AND PAIR BINDING 


Spin gap represents singlet to triplet excitation energy and is defined as: 


A,=E{N,Sj = l)-E{N,Sj = 0) 


( 1 ) 





(a) (b) (c) 

FIG. 1: (Color online) Spin gap extrapolation to thermodynamic limit at finite electronic filling n = 23/24. These plots are 
for different hopping modulation 5t and segment size Si = 2,4,6 in (a), (b) and (c) respectively with exchange coupling is 
J/t = 0.8 
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FIG. 2: (Color online) Binding energy extrapolation to thermodynamic limit at finite electronic filling n = 23/24. These plots 
are for different hopping modulation St and segment size Si =2,4,6 in (a), (b) and (c) respectively with exchange coupling is 
J/t = 0.8 


where E{N, S'J) is the ground state energy with quantum numbers N = nL and total z-component of spins S'J. Fig.[T] 
shows the finite size scaling of spin gap. We choose finite electronic density n = 23/24 and exchange coupling J/t = 0.8 
for different segment size Si =2,4 and 6 as indicated. In each case, the usual t—J model i.e., hopping modulation 
St = 0 gives zero spin gap. As we introduce finite St, the system develops finite spin gap and vanishes again at higher 
St. Similarly corresponding binding energy: A(, = E{N ± 2,S'J = 0) + E{N,S'^ = 0) — 2E{N ±1,S'^ = ±1/2) 
extrapolations are shown in Fig. [51 The spin gap and binding energies are found to be almost equal. 


LUTTINGER PARAMETER 


We find the transition from metallic to SC phase by calculating the Luttinger parameter Kp < 1 for metallic and 
Kp > 1 for superconducting phase. In thermodynamic limit Luttinger parameter is calculated from the slope of the 
structure factor for the density-density correlation at wave vector q 0 limit d, 13 • 


Kp 


= TT lim 

g—/0+ 


W{q) 

Q 


( 2 ) 


where for finite systems in numerical simulations q = 27r/L, L being the system size and the Fourier transform of 
correlation, 


W{q) = - {ni){nt+i))e 

The accuracy of Kp suffers in this approach due to the difficulties in accurate calculation of density-density correla¬ 
tions, specially at large distances and subsequently the Fourier transform at small q. So we calculate structure factor 
directly in momentum spaced extract Kp in thermodynamic limit by finite size scaling. 

For this we define an operator in momentum space: 

K{q) = ^(4'o|h((7)h(-(7)|4'o) (3) 

where n{q) = ^ ; Vc = {L + l)/2 being the middle position of the chain. Note that N{q) and W{q) 

1,(7 

are different, but becomes identical only in thermodynamic limit. So finally we get Kp by targeting not only the 
ground state 4>o but also the state h(( 7 )|'I'o) as: 

Kp = lim —N{2'k/L) 

L—¥00 2 


( 4 ) 
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FIG. 3: (Color online) The Luttinger parameter extrapolation to thermodynamic limit for a system with two holes i.e., n —^ 1 
limit at hopping modulation 5t = 0.2. These plots are for different exchange coupling J and segment length Si as indicated. 
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FIG. 4: (Color online) The Luttinger parameter extrapolation to thermodynamic limit for finite density n = 11/12. These 
plots are for different exchange coupling J and segment length Si as indicated. 


Fig. [3] shows the finite size extrapolation to thermodynamic limit for a system with two holes i.e., n —>■ 1 limit at 
hopping modulation 6t = 0.2. We see for each segment length Si the Luttinger parameter Kp in thermodynamic limit 
goes above Kp = 1 value where the metallic to superconducting transition happens. The critocal value of exchange 
coupling J'^ at this transition can be as small as ~ 0.35 as seen in S'/ = 6 system. The similar extrapolation of Kp 
at finite density n = 11/12 and for different S/ system are shown in Fig. S) In this way we have calculated Kp in the 
thermodynamic limit for different n and S/ to obtained the n — J phase diagrams shown in the main text. 


INSULATING LINE PHASES 

Due to the hopping modulation, the electrons tend to form some super structure in each segment causing insulating 
behavior at some commensurate fillings for different Si systems as shown in n — J phase diagrams in main text. Here 
we show charge gap Ac extrapolation to thermodynamic limit to show the melting of insulating phase to metallic 
state as we increase J. 

As we discussed in detail in the main text that for Si = 2 system at n = 1/2, the competition between t— and 
J—dimerization causes the melting to metallic states as we increase J. Here we discuss insulating lines for S'/ = 4 
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FIG. 5: (Color online) Typical charge gap extrapolation to thermodynamic limit for commensurate densities for Si = 4 system 
as we increase J. 


(Fig. 3(b) in main text) system in detail with the help of charge gap extrapolation. At n = 1/2 we notice that the 
Si = 4 system is insulating for all values of J up to phase separation. This is because two t—dimers are formed in 
each segment for small J and accommodate two electrons in bonding states of each t—dimer at this commensurate 
filling giving a Mott insulator. When we increase J, these two electrons already in bonding states lower the energy 
by pairing themselves. So this insulating phase cannot melt away with increasing J like the Si = 2 system at n = 1/2 
(discussed in main text). The insulating line phases are captured by calculating the charge gap defined as: 


A,= 


E{N + 2) + E{N - 2) - 2E{N) 


(5) 


Fig. Oja) shows the extrapolation of charge gap to thermodynamic limit at n = 1/2 for Si = 4 system. The charge 
gap remains finite for all values of J upto phase separation. However the melting of insulating states is possible at 
n = 1/4 and 3/4. At small J, one electron per segment (i.e., n = 1/4) occupies the lowest energy state formed 
by superposition of two bonding states already formed by t—dimerization. Also for n = 3/4, one electron sits on a 
t—dimer formed on middle two sites and other two electrons occupy the ends sites of each segment. As we increase J, 
the tendency of J—dimerization destroys both the insulating states at n = 1/4 and 3/4. Fig. [51(b) and (c) show that 
the charge gap vanishes as we increase J. It is intuitive that odd number of electrons in a segment is not suitable for 
effective pairings and lead to melting of insulating states as we increase J. Even number of electrons on the other 
hand can form an effective pairing which prevents the melting of insulating states as shown for n = 1/2 in S'; = 4 
system. 

Fig. 3(c) of main text depicts the insulating lines at different fillings for S; = 6 system. At n = 1/2, three electrons 
occupy bonding states in three <—dimers formed in each segment for small J leading to a Mott insulator. Again like 
Si = 4 system at n = 1/4 and 3/4, the odd number of electrons in each segment can’t make effective pair. So with 
increasing J we have melting of insulating state at n = 1/2 due to the competition of t— and J—dimerization. The 
similar argument holds for n = 1/6 and 5/6. But with even number of electrons at n = 2/6 and 4/6 form effective 
pairs as discussed in 5; = 4 system. So the system continues to be insulating upto phase separation. 

Finally in the phase separation the interaction between electrons becomes very large and the electrons forms 
antiferromagnetic domains such that the system separates into particle- and hole-rich region. We determine the phase 
separation boundary by calculating inverse compressibility defined as: 


X-' = ^{E{N + 2) + E{N - 2) - 2E{N)) 


( 6 ) 


After extrapolating energy to thermodynamic limit, we calculate x~^ as a function of J and obtain the phase separation 
boundary. Fig. |6|shows typical plot of x~^ as function of J at n = 1/4,1/2 and 3/4 for Si = 4 system. We consider 
the crossing of x~^ at zero line as onset of phase separation. 
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FIG. 6: (Color online) Inverse compressibility x 


as a function of J at different electron density for Si 


4 system 
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